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ABSTRACT 

The proximity prohle in the spectra of z ~ 3 quasars, where fluxes extend 
blueward of the He ll Ly a wavelength 304(1 + z) A, is one of the most important 
spectral features in the study of the intergalactic medium. Based on the HST 
spectra of 24 He II quasars, we hud that the majority of them display a proximity 
prohle, corresponding to an ionization radius as large as 20 Mpc in the source’s 
rest frame. In comparison with those in the H I spectra of the quasars ai z ^ 
the He II proximity effect is more prominent and is observed over a considerably 
longer period of reionization. The He II proximity zone sizes decrease at higher 
redshifts, particularly at ^ > 3.3. This trend is similar to that for H I, signaling 
an onset of He II reionization at z > 4. 

For quasar SDSS1253+6817 {z = 3.48), the He ll absorption trough displays a 
gradual decline and serves as a good case for modeling the He II reionization. To 
model such a broad prohle requires a quasar radiation held whose energy distri¬ 
bution between 4 and 1 Rydberg is considerably harder than normally assumed. 
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The UV continuum of this quasar is indeed exceptionally steep, and the He II 
ionization level in the quasar vicinity is higher than the average level in the in- 
tergalactic medium. These results are evidence that a very hard EUV continuum 
from this quasar produces a large ionized zone around it. 

Distinct exceptions are the two brightest He II quasars at z ~ 2.8, for which 
no significant proximity profile is present, probably implying that they are very 
young. 

Subject headings: Intergalactic Medium — Quasars: General — Ultraviolet: 
General 


INTRODUCTION 


The epoch at redshift ~ 3 marks the peak of qua sar formation (Richards et al.ll2006h 
as well as the reionization of the intergalactic helium flMeiksin I2009L M09 hereafter). As 
the reionization of singly ionized helium requires high-energy photons above 4 Rydberg, the 
hard UV background radiation field in the vast intergalactic space, commonly referred to 
metag alactic”, at z ~ 3 is believed to be ph otons originating from quasars instead of 


as 


hot stars (IHaardt fc Madadll996l: iMeiksin 1120051 1. Because helium is difficult to ionize and 
readily recombines, its opacity at this redshift range is considerably higher than that of 
hydrogen. 

Most of our knowledge of the intergalactic helium at these h igh redshifts is basec 


upon a few bright quasars with clear lines of sight : QS00302—003 flJakobsen et ah 


19961; 

lUechner et al. IboObh and HE2347—4 

o 

o 

Zhene et al. 

I 2 OO 4 ; 

Shull et al. 

2010) 


[l9^: 


York et ah 120001 1 and Galaxy Evolution Explorer flMartin et ah 1200511 . the number of 


known quasars with unobscured sightlines to their He II Ly a wavelengths has increased dra¬ 


matically from only three in mid-1990s to more than 50 fjSvphers et al.ll2009a : 


2009b 


Worseck et ah 1 1201111 . Of these quasars the highest redshift is 3.93 fjSvphers et al 


2009b|). The Gosmic Origins Spectrograph fGOS. lGreen et al. 1120121) instrument aboard the 


Svphers et ah 


HST has greatly enhanced our ability to probe multiple lines of sight and investigate the 
He II absorption features at higher signal-to-noise (S/N) ratios. In the hve years since the 
GOS installation, more than two dozen quasars have been conhrmed with the He II Ly a 
features. 


The ensemble of spectra of these quasars is now sufficiently large that we can begin to 
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overcome the cosmic variance and investigate the reionization proces s over a signihcant range 
of redshift. It is now well established that the Gunn-Peterson effec t flGnnn fc Petersonlll965l) 
in quasar spectra becomes signihcant only at ^ > 2.7 for He II flDavidsen, Kriss &: Zhene 


QQgI: Reimers et ah 1997Anderson et ah 1999 1 and at z > 5.5 for H I ( Becker et ah 


2001 


Fan. Carilli fc Keatinej 20061) . signaling the end of the respective reionization epochs. The 
intergalactic-medium (IGM) reionization process is lengthy, and little is known about the 
IGM evolution at much higher redshifts, i.e., z > 3.5 for helium, z >7 for hydrogen. 

The ionization properties of the IGM change signihcantly near luminous quasars, whose 
radiation enhances the ionization level in their vicinity. This “proximity effect” was hrst 
found in quasars at 2 ; ~ 3 where th 


spectr a declines in the quasar vicinity flMurdoch et ah Ill986 


1987b[l 


forest lines in 

their optical 

Garswell et ah 

1987; 

Tvtler 

in an expected residual hux 


blueward of the He II Ly a wavelength (jZheng fc Davidsen Il995l: iGiroux. Fardal &: Shull 
I 995 I: iMadau fc Reesll2000l). Such an absorption prohle has been o bserved in the H I spectra 
of most quasars at z > 5.7 (jWhite et ah II2003I : iGarilli et ah II 2 OIOI) . For helium, the presence 
of proximity prohles was ambiguous, as there were only a handful of known He II quasars, 
and some of them do not display a proximity prohle. 


The He II Ly a absorption troughs may be complex (IMadau &: ReesI I 2 OOOI) . as their 
shapes are affected by several factors: a proximity prohle blueward of the He II Ly a 
wavelength, a damped absorption prohle that extends to the wavelengths redward of it, and 
possibly another damped absorption prohle from an associated absorber. In this paper, we 
report our analysis of the proximity effect using a large sample of quasars with conhrmed 
He II Ly a features. The distances quoted i n this paper a re proper distances in the quasar’s 
rest frame, calculated using the codes of iHogg (jl999l) . with the standard cosmological 
parameter values: Gm = 0.3, Ha = 0.7 and h = Ffo/lOO kms“^ Mpc“^ = 0.70. 


2. OBSERVATIONS AND ANALYSES 

We investigate the HST spectra of 24 quasars in which the He II Ly a absorption feature 
is present at the anticipated wavelength and is not severely contaminated by geocoronal 
emission. Our database consists of three parts: (1) GOS spectra of four quasars at z > 3.4 
from our program (GO 12249: PI Zheng), see Tabled} (2) archival GOS spectra of seven 
quasars (GO 11742, 13013: PI Worseck); and (3) published spectra of 13 quasars, taken with 
GOS and STIS. Six other quasars are excluded because their He II break is contaminated 
by geocoronal Ly a or O I A1302 emission. Table [2] lists the eleven sources in parts 1 and 2; 
and Table E] lists the 13 sources in part 3. 
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Optical spectra of these quasars are useful as they provide information of the H I absorp¬ 
tion counterparts. We retrieved 15 SDSS spectra of the He II quasars, the VLT/UVES spec¬ 
tra of quasars PKS1935—692 and HE2347-4342, and the Keck/HIRES spectra of QS00302—003 
and HS1700-I-6416. In total, we have the optical spectra of 17 quasars, among which four 
are at high resolution. 


2.1. COS Spectra of Four Quasars 


We carried out HST/COS spectroscopy of four quasars: SDSS1253-I-6817, SDSS1319-I-5202, 
SDSS1711-I-6052, SDSS2346—0016 (GO program 12249). They were selected as a bright sam¬ 


ple (the UV continuum level at > 5 x 10 ergs ^ cm ^ A at z >3.4 with confirmed 


He II Ly a absorption brea ks with the ACS prism data (jZheng et al.l 120081: ISvphers et ah 


2009a : Svphers et ah 2009bh . The new COS data provide considerably higher S/N ratios 
and spectral resolution. The HST/COS observations of these four quasars were performed 
between 2010 November and 2011 December. The G140L spectra cover a wavelength range 
of approximately 1100-2150 A with a pixel scale of 0.08 A and resolution R ^ 2000. The 
source fluxes were extracted using a small slit width of 25 pixels , which is considerably 
smaller than the default size of 57 pixels (0".59, Svphers et ah 2012 : Svphers fc Shull 2013). 
To determine the level of geo coronal-line contamination and scattered light, we reduced 
portions of the orbital night separately and compared with the full set. Only in one case 


(SDS S2346—0016) the airglow is significant even during the orbital night (iFeldmau et ah 
19921 1. and the dataset taken in 2011 December has to be excluded. 


The spectrum of SDSS1319-I-5202 reveals a Lyman-limit system (LLS) at z ~ 0.7 (~ 
1550 A in the observed frame), and the source flux below 1550 A is considerably suppressed. 
Based on low-resolution prism data {R ~ 105 at 1460 A), we suggested an associ ated damped 


absor ption system whose effect extends redward of the He II Ly a wavelength fjZheng et ah 
20081) . With the improved wavelength accuracy of our higher-resolution COS spectrum. 


the absorption profile is conhrmed to straddle the He II Ly a wavelength and cannot be 
explained by a simple proximity profile. The COS spectra of these four quasars are shown 
in Figure [T] 

To study the absorption profiles, it is essential to accurately determine the quasar 
systemic redshifts. We re trieved the o ptical spectra of these quasars from the SDSS and used 
the IRAF task specfit flKriss Ill994j) to ht both the optical and UV spectra with multiple 
components. For the COS spectra, a power-law continuum, and a He II Ly a emission line 
were used to ht the data in wavelength redward of the He II Ly a break. For the optical 
spectra, we used a power-law continuum, multiple emission lines, including Ly a (narrow 
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+ broad), N V, O I, Si IV, C IV, He II fIVanden Berk et al. 1 120011) . and a set of absorption 
templates that simnlate the accnmnlated IGM absorption and the proximity profiles. The 
systemic redshifts were verihed nsing the O I A1304.22 emission line as low-ionization lines 
are believed to be less affected by the local motion of broad-line regions. For onr sample, 
the redshifts a re generally consisten^ t with those in the SDSS Data Release Seven Qnasars 
Catalog (DR7, Schneider et al. 2010 1. except for SDS S2346—0016. We ob tained a spectrnm 
of SDSS2346—0016 nsing the TripleSpec instrnment 1 Wilson et al. 2004 ) on the ARC (the 
Astrophysical Research Consortinm) 3.5-m telescope. The observations were made on 2011 
November 10 and 16, nnder an ambient temperatnre of 1-3 C, and the total integration time 
was 10 honrs. As shown in Fignre[2l the [O III] A5007 emission is not visible. The Mg II and 
H/5 emission lines, while nnnsnally weak, yield a redshift consistent with the O I line. The 
COS spectra are of snfficient spectral resolntion to conhrm that that onr measnrements of 
the proximity profile are not affected noticeably by the small difference in qnasar redshifts. 


2.2. Archival COS Spectra of Seven Quasars 


We retrieved the archival COS/C140L spectra of 15 other qnasars (CO 11742 and 
13013) via the Miknlski Arch ive for Space Telescopes and processed with the standard 
pipeline calcos flHodge Il201ll. v2.19). The difference between the pipeline extraction and 
onr narrow-window extraction (see 112 .ip is mainly in the residnal flnx level at wavelengths 
far away from the qnasar, bnt this difference has no effect to onr stndy. For QSO0233—0149, 
the geocoronal O I A1302 line is close to the He II Ly a wavelength, therefore we nsed the 
orbital-night portion of the data. We also nsed the night portion of the QS00916-I-2405 
data as geocoronal O I A1356 emission is near its He II proximity prohle. Fonr of these 
qnasars were exclnded becanse of severe contamination by geocoronal emission. Fonr more 
were exclnded as the flnx breaks in their spectra are observed at wavelengths considerably 
longer than that for He II Ly a. Most likely, these are attribnted to partial He II LLSs 
along the lines of sight. The details of the seven archival spectra nsed are listed in the lower 
parts of Tables d] and dl below the details for onr prime COS spectra listed at the top of 
the tables. Fignre [3]plots the spectra of seven qnasars that we inclnde in onr sample, with 
contaminating lines flagged with the Earth symbols. 


2.3. Archival Data in the Literature 


A number of the COS spectra of He I I quasars have been published fIShull et al.l 12010 


Svphers et al.l 120121: ISvphers fc Shnlll 120131) . most of which were observed with the G140L 
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grating. We estimated the proximity sizes in the literature. For most of them, their systemic 
redshifts are from the SDSS catalog. The references for their redshifts are listed in Tables 3 
and 4. One STIS spectrum of quasar PKS1935—692 is used. 


2.4. Measurement of Proximity Zone Size 


The study of He II proximity prohles is more challenging than that of H I at 2 ; ~ 6 as 
they are more complex. We measured a proximity zone from the He II Ly a wavelength to a 
point where the flux drops to below propagated errors as calculated in a bin of approximately 
1.5 A. We used these nominal bin sizes and combined the redshift errors to estimate the 
uncertainties in zone sizes. A difference of 1.5 A corresponds to a redshift difference of 
0.005 for He II Ly a. This definition of proximity zones is based on an assumption of zero 
continuum flux beyond a proxirnity zo ne, as implied by the high He II opacity at z > 2.8 
f Shull et ah 2010 : Svphers et ah 2011a). 

An exception is the quasar HS1700+6416 {z = 2.751), where the con tinuum flux level 
blueward of the He II Ly a wavelength is known to be nonzero. In Figure 8 of ISvphers fc Shull 
(120131) . the He II Ly a absorption region is plotted along with a Keck/HIRES spectrum. We 
estimated the flux level in two regions where no strong H I absorption is present: 2.652 < 
< 2.653 and 2.6545 < 2 ; < 2.655. The average He II Ly a opacity there was estimated 
as r = 0.95. E xtrapolating this value to ^ = 2.73 using a theoretical scaling relation 
Teff oc (1 + z)^'^ flFardal. Giroux fc Shull Ill998l) . the normalized flux level there should be 
0.36 (r = 1.03). This is approximately the flux level around z ~ 2.728 and 2.739, where 
no strong H I absorption features are present, and a potential flux excess near these two 
redshifts is marginal at best. It is likely that the proximity zone in HS1700+6416 ends at the 
absorption feature of 2 ; ~ 2.744. A generous estimate would mark the zone down to 2 ; ~ 2.73, 
where there is no signihcant H I absorption line. We therefore list the zone size as 4±4 Mpc 
in Table [3l to reflect the two possible values of 2 and 6 Mpc and other uncertainties. 


All broad proximity prohles are embedded with strong absorption features. At a rela¬ 
tively high hux, these discrete absorption lines can be easily identihed. However, at a low-hux 
end, they may cut off the residual hux and cause an underestimate of the proximity zone 
size. When an optical spectrum is available, we checked whether a strong H I absorption line 
is present near the wavelength that corresponds to the He II proximity zone’s endpoint. In 
six cases in Tables [2] and [3], no such a H I absorber is identihed within 6 A of the endpoint’s 
wavelength in optical spectra {dz/{l + z) = 0.001). In the other 11 cases such H I absorption 
features are found near the proximity-zone endpoint, and for the other seven quasars no 
optical spectra are available. Therefore the possibility is real that their proximity zones may 
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be underestimated. In these 18 cases we increased the measurement errors by adding a term 
of 3 Mpc {dz/{l + z) = 0.0035 at ^ = 3.5), approximately the nominal width of a strong 
absorption feature in G140L spectra. The quasar I1E2347—4342 does not show a proximity 
zone, and there is no absorption feature in the UVES spectrum near the systemic redshift 
that can account for a potential He II LLS. 


2.5. Correlation with Redshift and Lnminosity 


Given the sample size, it would be informative to study the evolution of He II proximity 
zones with redshift and luminosity, which has been reported for H I at 2 ; ~ 6 f 1 13.81) . We plot 
the distribution of proximity-zone sizes over redshift and luminosity, respectively, in Figures 
[5]and[6l The redshift scales in the top and bottom panels of Figure [5] are set to reflect the 
same ratio of cosmic time. Our results demonstrate that the proximity effect is common 
among He II quasars, and their proximity-zone sizes can be considerably larger than those 
in the H I quasars at z ~ 6. The most significant trend seems to be at z > 3.3, where the 
zone size decreases towards higher redshift. This trend may be understood in terms of an 
increasing He II fraction and signals the onset of the intergalactic He II reionization. 


The two most luminous quasars, HS1700-I-6416 and HE2347—4342, are located at lower 
left in the top panel of Figure[H]and clearly d isconnected from the rest of the sample. Recently 
a bright quasar at z = 6.3 was discovered (jWu et ah 1120151) : with a similar luminosity and 
a moderate zone size, this quasar is also at odds with the model prediction that luminous 
quasars display a larger proximity zone. 

We ran statistical tests to check potential correlations. The tool is t he IRAF astronomi- _ 

cal su rvival analysis package, which allows upper or lower limits as inputs (jlsobe. Feigelson fc Nelson" 
19861) . A moderate anti-correlation between zone sizes and quasar redshifts is confirmed: 
the Spearman correlation coefficient p = —0.55, rejecting a no-correlation hypothesis at 99% 
level. Proximity-zone sizes do not appear to be correlated with quasar luminosities: the 
Spearman correlation coefficient between them is 0.09, suggesting a 65% probability for no 
correlation. Even without the two data points at the high-luminosity end, the correlation is 
still poor. 










3. DISCUSSION 


3.1. Uncertainties in Zone Size 


A major source of uncertainties in our measurements of proximity zones is the sys¬ 
temic redshifts of q uasars. The SPSS qu asars have redshift uncertainties on the order of 
dzl{V -f 2 ;) ~ 0.002 flHewett fc Wildll2010l) . which are derived mainly at z < 0.8 when nar¬ 
row emission lines are present. We verihed the redshifts of SDSS quasars in our sample 
with O I A1304 emission in the optical spectra. O I A1304 emission is a weak feature and 
consists of mu l tiple components , which carry a statistically weighted average of 1303.49 A 
fjMortonlll99ll: IShull et ahl 120101). Furthe rmore, O I A1304 emission may be blended with 
Si II A1305.42 (IVanden Berk et ah 1 1200 if) . Assuming a solar abundance, the weighted av¬ 
erage for O U-Si II is 1304.22 A. We adopted this value with an uncertainty of 1 A, which 
is comparable to the statistical errors of 5 A (in the observer’s frame) in our htting to the 
SDSS spectra. Combining these terms and converting into the observer’s frame for He II 
Ly a, we estimated uncertainties of approximately 2 A in the UV spectra of our sample. 
This term corresponds to dz/{l + z) = 0.0015, 450 kms“^ in velocity dispersion, and 1.2 
Mpc in proper distance at z = 3.5. We verihed the redshifts of SDSS quasars using hts to O I 
A1304 emission when possible. In a number of cases, O I A1304 emission is too weak to be 
useful in redshift estimates, we then used the redshift values from DR7. For the non-SDSS 
quasars in Tables [2] and [3l we assumed a redshift error of 
unless quoted explicitly. 


800 kms in velocity space. 


Another source of uncertainties in our estimates of the proximity-zone sizes arises from 
their ending points in the low-hux regions. Our data are mainly moderate-resolution spectra 
with limited S/N levels. For the COS G140L data, this nominal uncertainty is 1.5 A, which 
corresponds to approximately 0.9 Mpc in proper distance at z = 3.5 and 1.1 Mpc at z = 3.2. 
As discussed in 1 )2.41 the frequent occurrence of strong absorption features (H I column 
density ~ 10^^ cm“^) adds a term of uncertainty in many cases. 


Other rare effects, such as an infalling absorber could affect t he redshifts, but a re unlikely 
to be thick enough to account for the lack of a proximity effect (jShull et al.ll2010h . 


3.2. Effect of Intervening Hydrogen Lyman Limit Systems 

Quasar spectra display the signature of numerous IGM components, and density fluc¬ 
tuations have a signihcant impact on the He II ionization process. LLSs of column density 
> 10^^ cm“^, both in H I and low redshift and He II at high, produce random flux cutoffs. 
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At 2 ; ~ 0.5, a H I LLS may mimic a He II absorption edge. Th e probability of encounter¬ 
ing s uch a system per unit redshift is dn/dz = 0.28(1 + (Riba udo. Lehner & Ho\™ 

2011 ), which yields a rate th at is consistent with other results f Storrie-Lombardi et al.fllQlM 


Stengler-Larrea et ah Ill995l) . Over a wavelength range of 15 A 10 Mpc in proper dis¬ 


tance) near He II Ly a, there is less than 2% probability of intercepting a H I LLS of optical 
depth r > 1. This value would be even smaller if we consider this effect only in the region 
near the proximity-zone’s endpoint. 


3.3. Effect of intervening He II Lyman Limit Systems 


For He II, the IGM components have considerably higher column densities and a filtering 
power. If the ratio of column densities rj = NheIi/Nhi ^ 100, even a moderate H I absorber 
may become a potential He II LLS greatly reduce the number of He Il-ionizing photons 
from the quasar, reducing both the growth of a He II ionization zone and the subsequent 
ionization rate. 


A simple estimate flMcQuinn et ah II2009I) suggests that a system becomes opaque (r > 
1) to the ionizing photons at 4 Rydberg if its H I col umn density > 1.3 x lO^^fhO/r?) 


cm 


-2 


Based on the high-resolution data at z zz 2.8 flKriss et al.ll 200 ll: IZheng et al.l 12004 


Fechneretah 


20061 ) . the mean value of rj is around 80, which is derived from a majority 
of weak components (column density Nht < 3 x 10( ^ cm“^). At 3 < z < 3.2, the mean 
Tj value may be around 200 fjSvphers fc Shn11 l2014l) . However, very high rj values are 


mostly associated weak lines and very difficult to confirm, particularly at z > 3 when the 
He II opacity is high . Evidence suggests that strong components display lower rj values: 
Fechner et ah ( 120061 ) found 77 A 15 at N^i > 3 x 10 ^^ cm ~^, which are supp orted by the 
data in another bright quasar (77 A 10 at Nhi > 10 ^^ cm~^. IZheng et ahl 120041 ) . Adopting a 
hducial value of 77 = 15 , only the systems of H I column density > 4 x 10 ^ ^ cm“^ may hlter 


out the quasar He II-ionizing radiation. From Table 3 of iKim et al.l (120131) the distribution 
can be expressed as d‘^n/[dNuidX) = 10® ®"^ The probability of encountering H I 

systems with Nhi > 4 x 10^® cm“^ within 1.5 A to the endpoint of a proximity zone (see 
^2.4p is estimated as ~ 0.02. 


Theoretical models, however, predict high 77 around 100, based on our knowledge of 
the quas ar’s EUV (extreme UV) co ntinuum, and a re-analysis of the average 77 values at 
X ~ 2.5 (iMcQuinn fc Worseck 1120141) does not hud evidence for the extremely low 77 values 
noted in previous studies. This recent analysis is based on G140L data, while the previous 
studies used high-resolution FUSE data. To check the realit y of low 77 at x > 3, we used 
the published results of QS00302—003 at a higher redshift. Svphers fc Shull ( 2014 ) used 
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a high-resolution Keck spectrum and COS/G130M spectrum to derive r] values with a bin 


In the range of 3.03 < z < 3.235 there are 21 H I absorption lines at 
flKim et ah 2002). At 2.76 < 2 ; < 2.93 and 3.25 < 2 ; < 3.267 there are 


size of dz = 0.00 
N > 10^^ cm-2 

additional eight absorption lines with re stframe equivalent widths greater than 0.3 A. From 
Figures 10-13 of ISvphers Sz Shulll (120131 ) we estimated 29 rj values: twelve of them are lower 
limits and 17 are measurements. For the twelve lower limits, their mean value is r; = 12. For 
the 17 measured r] values, their mean value is 13. 


If He II LLS are common, no broad proximity profiles should have been observed. We 
analyzed the SDSS spectra of our sample and identihed 40 strong H I absorption features 
in the proximity zones: nearly all see a significant flux at their blueward wavelengths. If 
the r] values are around 100 for these absorbers, they should display a sharp flux cutoff. 
Furthermore, in the SDSS spectra of hve quasars, no strong H I absorbers are found near 
the wavelengths that correspond to a proximity-zone endpoint (Tables |2] and [3]). The upper 
limit to the H I column density would be approximately 10^^ cm“^, assuming a minimum 
equivalent width of 0.3 A in the restframe and a nominal velocity dispersion of 30 kms~^. If 
an endpoint of the proximity profiles in these quasars is attributed to an assumed He II LLS, 
the Tj value may be 8000 or higher, which is not supported by observations. We therefore 
conclude that, while such a possibility cannot be ruled out, it is unlikely that proximity 
prohles are signihcantly altered by dense He II LLS. 


3.4. He II Ionized Zone 



a quasar, and the quasar radiation is isotropic. The total quasar ionizing luminosity in the 
restframe, in units of photon rate, is calculated as = 47rD|(l -|- z)~^ where 

Dl is the bolometric luminosity distance, fi, is the flux density in the observer’s frame and 
z/q is the observed frequency corresponding to 4 Rydberg. The Stromgren radius is then: 


Rs = 


3W 


1/3 


AnarienHer 


73 


N ^/10 


57^-1 


1/3 


A2 


1 + ^ 


-2 


Mpc 


( 1 ) 


where a is the Case-A recombination coefficient of He III ions, Ue the electron density, 
and A the local baryon overdensity in units of the mean IGM density, and i?s the proper 
distance. The mean density of 1.85 x 10“^ cm“^ at 2 ; = 3.5 is derived from its value at ^ = 0 
( Komatsu et ah 2011 : Svphers 2010). 
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The value given in Equation (1) is an upper limit; real ionized zones are considerably 
smaller for at least two reasons. First, this hypothetical limiting radius grows faster than the 


ioniz a tion front expands due to the decreasing IGM density as the Universe expands flShapiro 
19861: iDonahue fc Shull 1 119871) . Secondly, for intergalactic helium, the recombination time 
scale is so long, (cme)“^ ~ 1 Gyr, that it exceeds the nominal quasar ages. As a result, 
the IGM around quasars would not reach a balance to form a Stromgren zone. If the gas 
is denser than average (A > 1) or clumpy, and if the quasar age is old, the recombination 
rate may become comparable to the ionization rate. A proximity zone does not mean that a 
Stromgren sphere has been established, in the sense of numbers of recombinations balancing 
numbers of ionizations. This calculation could just indicate that the gas within the sphere 
is over-ionized by the presence of the quasar radiation compared with that in the general 
IGM. The timing criterion for the ionization fraction to stay in ionization balance is only 
that the age be longer than (ionization rate)“^ for highly ionized gas (M09). 

Assuming an average IGM density, the opacity is estimated t o be on the order of 


2 X 10^ [(1 -f- for hydrogen and 4 x 10^ [(1 -|- for He II l MiraIda-Escudall998 


Svphers 11201011 . ISvphers et al.l (l2011bl ) provide a more accurate formula that can be written 
for the chosen cosmological parameters: 


T'GP.Hell — 3.49 X 10^ A a^Hell 


Yp 


0.2486 


1 + ^ 


3/2 


1.0209 


_1 + 0.0209- (4/(1 +;2))3_ 


( 2 ) 


where XHeii is the He II fraction and Ij, is the helium mass fraction. 


In the vicinity of a luminous quasar, the ionizing radiation held is enhanced over 
any metagalactic ionizing backgronnd radiation, producing a cosmic “bubble” that may 


be observed in t 


re quasar spectrum as a proximity effe ct. This effect s hould apply to H I 


flMurdoch et al. 


1995 


1986; 


Baitlik, Dunc an fc Ostrikeiill988h as well as He II fjZheng fc Davidsen 


Hogan. Anderson fc Rugersil997 ) absorption. However, the overall He II optical depth 


could be as high as ~ 4 x 10^ if the intergalactic helium is in a form of singly ionized state 
and remains poorly constrained at z > 3.3. 


In the wavelength region 1325-1340 k{z = 3.362—3.411) of the spectrum of SDSS1253+6817, 
the effective optical depth is Tueii = 4.101q; 33 9^^ conhdence level. If we used a wider slit 

width as the calcos pipeline default, r = 5.10lo;52 at 68% conhdence level, or r = d.lO/^gy 
at 95% level. It is therefore possible that the hux extends beyond 12 Mpc from the quasar, 
suggesting that the metagalactic He II-ionizing radiation held is extremely weak but present, 
or the quasar has a He Il-ionizing luminosity considerably stronger than the extrapolation 
from a nominal power law, or the ionization level has not yet recovered after an earlier 
episode of an enhanced quasar luminosity. 
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3.5. Effect of Quasar Age 

The proximity profile is subject to the quasar age. A large zone size requires sufficient 
ionizing photons from the quasar, even after a long period, if the surrounding intergalactic 
helium was not fully ionized, and a large number of photons are needed to ionize it. For 
a young quasar, the proximity prohle would display a sharp cutoff. Several factors affect 
the proximity prohle: (1) the IGM density and its huctuations; (2) the quasar luminosity; 
(3) the quasar age; and (4) the recombination process. For the case when the quasar age is 
considerably smaller then the recombination time scale, the total number of He 11-ionizing 
photons over a lifetime tg is equal to that of the ionized helium atoms in the quasar vicinity: 


-^tQ 


yATTHHeXHell ) 
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'( iv^/io5V^)(tQ/ioV )' 

3^HeII A 


1/3 


1 + ^ 


-1 


Mpc. (3) 


We used the proximity prohle in SDSS1253+6817 to model the evolution of intergalactic 
helium. While other cases of signihcant proximity ehect exist, the prohles are often shelf¬ 
like, possibly due to the duty cycle of the quasar or the transverse proximity ehect by 
other quasars (ISyphers fc Shull 1120141 ). The G140L spectru m, after a reddenin g correction 
of E{B — V)= 0.02, Ry = 3.1, and the extinction curve of iFitzpatrick fjlOOOl) . displays a 
threefold increase in hux from 1950 to 1400 A, which is steeper than that in QS00302—003. 
Since no LLS break is obvious in this wavelength range, we htted the dereddened continuum 
between 1400 and 1900 A with a simple power law and a potential partial LLS between 2000 
and 3000 A. The best ht suggests j3 = —2.09 (/a oc A^) and a LLS of r = 0.83 at ~ 2450 
A. If no LLS is assumed, the best-ht power law has [3 = —2.85. Both htting results suggest 
a similar continuum level of (4.6 ± 0.5) x 10“^® ergs“^cm“^A at the He II Lyman-limit 
wavelength of 1020 A. The estimated total photon rate at 4 Rydberg is (5.9±0.6) x 10^® sec“^. 
As shown in the next section, this hard continuum is not strong enough to produce the 
observed broad proximity prohle if the intergalactic helium is not fully ionized. 


3.6. Time Dependence of Ionization 


The ionized zones represented by Equations (1) and (3) are extreme cases; the real 
ionization by a quasar is a gradual and slower process, which can be described by aii 
analytical expression d erived from a simplihed differential equation flGen fc HaimanI 12000 : 
Svphers fc Shull II2014I) for a hxed IGM temperature. A more accurate evaluation of the ion¬ 


ization fraction requires simultaneously solving for the gas temperature and the ionization. 
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We computed the time dependence of He II ionization aronnd a qnasar, as ontlined in M09. 
Specifically, we solved the time-dependent spherically symmetric set of conpled ionization 
rate eqnations for nniformly distribnted hydrogen and helinm (Eqnations 41 and 46 of M09) 
on a radial grid, inclnding the photoelectric attenuation by H I, He l and He II of the radi¬ 
ation from the quasar. We solved these equations simultaneously with the photoionization 
heating and cooling equations as described in M09 (sec.HI.B) for a pure hydrogen/helium gas 
of cosmic abundances. This calculation used the atomic rat es in M09, except for ado pting 
the electron excitation and ionization cooling rate of H I from IScholz fc Walters (1199111 . Our 
calculations assumed Case A recombination and included Compton cooling off the cosmic 
microwave background and adiabatic expansion losses in the gas. We considered a sequence 
of quasar activation redshifts corresponding to He Ill-region ages ai z = 3.48 of 10, 20, 30, 50, 
and 100 Myr. To match the observed proximity profile, the He Il-ionizing continuum must 
be very strong: = 6 x lO^^erg s“^Hz“^ at 228 A in the restframe, which is approximately 

three times of the estimated level from observations. It would imply an EUV continuum 
that is even steeper than QS00302—003 fjSvphers fc Shull 11201411 . Lower-luminosity models, 
such as a He Il-ionizing continuum extra polated from the flux at He II Ly a w avelength 


with a nominal power-law index of —1.7 flZheng et ah 119971: iTelfer et al.l 1200211 were not 


able to match the size of the ionized region and the flux simultaneously. We emphasize 
that the measured absorption signal is co-temporaneous with the quasar observed. Since 
the characteristic time to reionize the gas and maintain ionization equilibrium is on the 
order of millions to tens of millions of years, depending on the distance from the quasar, the 
luminosity of the quasar that produced the proximity zone may have been larger than the 
observed value; we have few constraints on quasar variability on million year timescales. 

We estimated the black-hole mass as 2.4 x 10® Mq, using the C IV line width, the un¬ 
derlyi ng continuum flux in th e SPSS spectrum and the formula in IVestergaard fc Peterson 
(120061) . Following the work of ISteinhardt fc ElvisI (1201 ll) . we derived the rest-frame luminos¬ 


ity near the Lyman limit as 2 x 10^^ ergs“^, which is below the Eddington limit of 3 x lO'^^ 
ergs“i. A He Il-ionizing metagalactic background was not assumed. In our simulations 
of the ionization zone, a hydrogen-ionizing and He l-ionizing metagalactic background was 
turned on at ^ = 7. Figure [7] shows the He II ionization level, the IGM temperature and the 
normalized observed He II flux at different epochs after the quasar’s birth. 


Shortly after the quasar’s birth, a sharp ionization front expands over time, as described 
by Equation (3). At each radius and each time interval, the temperature is computed from 
the energy equation and used to evaluate the atomic rates. Helium within the sphere is 
highly ionized, as the ionizing flux is strong and recombination is not effective. At later 
epochs, when the ionization front reaches large distances, ionization balance between the 
geometrically-diluted quasar ionization flux and recombinations produces a more gradual 
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rise in the He II fraction with distance. As shown in Figure [TJ a smooth decline of flux may 
be observed after ps 30 Myr. 

The high luminosity assumed in this simulation and that in Equation (3) probably rep¬ 
resents an upper limit as high-energy photons from the quasar need to ionize a large volume 
of singly ionized helium. This requirement would be eased if a He Il-ionizing background 
held existed prior to z ^ 3.5. As Equation (j3]) shows, the ionizing luminosity may be re¬ 
duced by a factor of three to 6 x 10®® sec“^ if the He II fraction XHeii is 30%. In either 
ca se, the ionizat i on fro nt is likely moving nearly at the speed of light, as shown in Figure 


[3 IWhite et ah (120031) studied the time-retardation effects on the observed proximity-zone 


size. Adapting their equations to He II, the speed-of-light corrections are important to the 
observed size until times large compared with ~ 100 Myr for an assumed boosted quasar 
He II ionizing photon luminosity of 1.9 x 10®^sec“^, where tc = (3iV^/47rnHeiiic®)^/^. The 
observed ionization zone sizes then increase slowly with time, as After an expansion 
time of 50 Myr, the ionization front in Figure [7] is restricted by the speed of light, but the 
ionization front velocity slows to 0.5c by 100 Myr. According to Equation (3), the age of 
the observed quasar producing the ionization zone after an expansion time of 50 Myr is 
tq = 12 Myr, which provides a minimal age to the quasar. 

If the IGM were highly clumped, the He II region would expand more slowly, as recom¬ 
binations would slow it down. The high quasar luminosity used in Figure [7] is required to 
match both the position of the He III ionization front and the He II Ly a flux level behind 
the front assuming a homogeneous IGM. If instead the IGM opacity is dominated by line¬ 
blanketing by the Ly a forest, the quasar luminosity need not be so high to match the flux 
level while still matching the position of the He III front. This would, however, require a 
somewhat older quasar. The model in Figure [7] demonstrates the extreme assumptions that 
must be made for the quasar spectrum to match both the flux level and the position of the 
He III front assuming a homogeneous IGM. Line-blanketing by the H e II Ly a forest is a 
more plausible explanation for the flux level flMadau &: Meiksin 1119941) . 


3.7. Helium Proximity Effect in the Literature 

In principle, a proximity profile should be present in every He II quasar spectrum, as a 
cosmic “bubble” would exist even without a metagalactic radiation field. At the beginning, 
ionizing photons do not travel very far, as the mean free path {aHeiinHeii)~^ is less than 
0.1 Mpc, where aneii = 1-6 x 10“^® cm^ is the photoionization cross section at 4 Rydberg 
and riHeii — 2.5 x 10“® cm“® the number density of singly ionized helium at ^ = 4, is 
less than 0.1 Mpc. High-energy photons from the quasars, however, penetrate considerably 
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deepe r into the IGM, raising its temperature flTittlev &: Meiksin 120071: iMcQninn et al 
20091) . These ionized spheres expanded and eventually overlapped, gradually completing the 
reionization of intergalactic helium. While it is anticipated that luminous quasars exhibit 
a broader proximity prohle (Equations (1) and (3)), it is not the case for the two brightest 
He II quasars at z ~ 2.8. The UV sp ectroscopic properties of HS1700+6416 (z = 2.751 


V=16.2) have been ext ensively studied flDavidsen. Kriss fc Zhendll996l: iFechner et al. 112006 


Svphers fc Shull! l2013l) . Its proximity profile is insignificant, extending at n iost 6 Mpc 


see 


1997; 

Kriss et ah 

2001 ; 

Shull et ah 

2004; 

Zheng et ah 

2004 

) has no proximitv effect observed. 


The re dshifts from O I A1304 and O III A5007 agree flReimers et al. 119971: ISvphers et al 
2011ar). There may be a strong infall i ng absorption syst em observed in the UV spectrum 


flFechner et al.ll2004 IShull et al.ll2010[l . IShull et al.l (120101) studied this infalling system and 
concluded that no satisfactory explanation exists for the absence of a proximity effect around 
this luminous quasar other than the possibility that it has only recently turned on (within 
the last Myr) 


Hogan. Anderson fc Rueers (119971) and iHeap et ahl (120001) discovered a He II proximity 
zone in the quasar QS00302—003 at z ~ 3.286 (V=17.5). Blueward of the break, a shelf 
of residual flux co rresponding to an opac ity r ~ 0.8 extends out about 15 A. The recent 


COS observations (jSvphers fc Shull 1120141) provide details in the proximity zone. The broad 
proximity prohle, extending up to 15 Mpc, may be attributed to a line-of-sight effect from the 
quasar itself, or to a transverse effect by another quasar Q0301—00 (z = 3.232) near the line 
of sight. In the latter case, the age of the ionization zone around Q0301—00 should be at least 
34 Myr. It is interesting to note that the He Il-ionizing continuum is considerably harder 
(power-law index —0.8) than previous work suggested: broad proximity proh les appear to be 


relate d to a hard quasar continuum. Similarly, in PKS1935—692 (z ~ 3.18, lAnderson et al. 


19991 ). a proximity shelf of hux extends blueward of the He II break by at least 20 A. 
This absorption feature displays a strong recovery void at 1246.5 A, likely produced by the 
radiation of a foreground source. 


3.8. Comparison with Hydrogen Proximity Effect 


The reionization of intergalactic h ydrogen ended by redshift ^ p s 6 , as characterized by 
a sharp increase in opacity at 2 ; > 5.7 (iFan. Carilli fc Keating] 120061 ). This trend is similar 
to what is observed for helium at z ~ 3, and it is instructive to compare these two major 
cosmic processes. The number of H I-ionizing photons of a quasar is considerably higher 
than that of He Il-ionizing photons. On the other hand, the IGM at z ~ 6 is denser than 
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that at z ^ 3.5. The recombination time scale of hydrogen is lower than helium. As a result 
of these factors, our estimate based on Equation (1) suggests that the proximity prohle of 
intergalactic hydrogen at 2 ; ~ 6 is somewhat weaker than that for He II at ^ ~ 3. 


Carilli et ah 


(j2010[) reported proximity prohles of hydrogen in the spectra of 27 quasars; 
their “near zones” (proximity) lie in the range of ~ 5 — 10 Mpc. They found a significant trend 
for a decrease in the near-zone size with increasing redshift, as evidence for the evolution 
of an increasing neutral fraction of intergalactic hydrogen toward higher redshifts. They 
also found that the near-zone sizes increase with the quasar UV luminosity, as expected for 
photoionization dominated by quasar radiation. We add eight more data points from recen t 


literature. One quasar is at the highest known redshift 2 : ~ 7.08 flMortlock et ah N 2 OIII) 


and another is luminous (M = —29.3, IWu et ah II2015I) . These two data po i nts ca rry a 


significant weight in the correlation analysis. The paper of IVenemans et ah (120131) does 


not list the near-zone sizes, and we made estimates from their Figure 4. A Spearman test 
over these 35 data points finds a correlation coefficient p = —0.59 between zone sizes and 
redshifts, rejecting a no-correlation hypothesis at 99.95% level. Another test between zone 
sizes and absolute magnitude finds p = —0.47, rejecting a no-correlation hypothesis at 99.4% 

level. 


The luminous quasar JOlOO-l-2802 displays a proximity zone of 8 Mpc. While this is a 
signihcant size, it does not scale well to luminosity (see the lower panel of Figure | 6 ]). Its 
normalized zone size, Rnorm = R at the restframe wavelength of 1450 A, is 

4 Mpc, below the average at the absolute magnitude M 1450 = —27. This result is actually 
consistent with the two luminous He II quasars that they do not show a signihcant proximity 
prohle. Maybe some other factors, such as the density and ionization level in the quasar’s 
vicinity or a young quasar age, play an important role. 


4. CONCLUSION 

We utilize a sample of 24 He II quasars to investigate the interplay between these 
quasars and the surrounding IGM. The ionization zone around He II quasars is often more 
prominent than that for H I at ~ 6. The large redshift range for these quasars allows us to 
gain insight into the IGM reionization over a long epoch before intergalactic helium became 
fully ionized. The proximity-zone sizes decline signihcantly at z > 3.3, and it is likely that 
helium reionization started well before 2 ; = 3.8. 

In the quasar SDSSI253-I-68I7, the source hux extends considerably blueward of the 
He II Ly a wavelength, suggesting a quasar age of > 12 Myr. The UV hux rises dramatically 














17 


from 1950 to 1400 A, suggesting an exceptionally hard EUV continuum. The t] value in the 
proximity zone is lower, consistent with such strong He Il-ionizing radiation from the quasar 
that produces a broad proximity zone. 

The two brightest quasars do not display a significant He II proximity profile. While 
this is at odds with model expectations, we notice that a luminous quasar does not display 
the largest H I proximity profile either. It is possible that these hyperluminous quasars are 
young, or they are surrounded by a dense IGM. 

It should be stressed that the observed proximity-zone size is not a direct measurement 
of the quasar lifetime, as the structure of the zone is affected by many possible factors, 
including the IGM density fluctuations and quasar variability over IGM ionization timescales. 
A quasar with small proximity zone may be considerably younger compared with the light- 
travel time across its proximity zone allowing for retardation effects. Further understanding 
of the quasar lifetime awaits improved IGM simulations that take these factors into account. 


We thank the anonymous referee for many thoughtful and instructive suggestions. 

Support for this research was provided by NASA through grant GO-12249 from the 
Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 

This work is based on observations made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555 and 
observations obtained with the Apache Point Observatory 3.5-meter telescope, which is 
owned and operated by the Astrophysical Research Gonsortium. 

Funding for the SDSS and SDSS-H has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Science Foundation, the U.S. Department of 
Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, 
the Max Planck Society, and the Higher Education Funding Gouncil for England. The SDSS 
Web Site is http://www.sdss.org/, 

This research has made use of the Keck Observatory Archive (KOA), which is operated 
by the W. M. Keck Observatory and the NASA Exoplanet Science Institute (NExScI), and 
the data obtained from the ESO Science Archive Facility under request number 154374. 



- 18 - 

REFERENCES 

Anderson, S. F., Hogan, C. J., Williams, B. F. & Carswell, R. F. 1999, AJ, 117, 56 
Bajtlik, S., Duncan, R. C. & Ostriker, J. P. 1988, ApJ, 327, 570 

Becker, G. D., Hewett, P. C., Worseck, G., & Prochaska, J. X. 2013, MNRAS, 430, 2067 

Becker, R. H., Fan, X., White, R. L. et al. 2001, AJ, 122, 2850 

Carilli, C. L., Wang, R., Fan, X. et al. 2010, ApJ, 714, 834 

Carswell, R. F., Webb, J. K., Baldwin, J. A. & Atwood, B. 1987, ApJ, 319, 709 

Cen, R. & Haiman, Z. 2000, ApJ, 542, 75 

Dall’Aglio, A., Wisotzki, L., & Worseck, G. 2008a, A&A, 480, 359 
_, 2008b, A&A, 491, 465 

Davidsen, A. F., Kriss, G. A. & Zheng, W. 1996, Nature 380, 47 

Donahue, M. & Shull, J. M. 1987, ApJ, 323, L13 

Fan, X., Carilli, C. L., & Keating, B., 2006, ARAA, 44, 415 

Fardal, M. A., Giroux, M. L., & Shull, J. M., 1998, AJ, 115, 2206 

Fechner, C., Baade, R., & Reimers, D. 2004, A&A, 418, 857 

Fechner, C., Reimers, D., Kriss, G, A. et al. 2006, A&A, 455, 91 

Feldman, P. D., Davidsen, A. F., Blair, W. P. et al. 1992, Geophy. Rev. Lett., 19, 453 

Fitzpatrick, E. L. 1999, PASP, 111, 63 

Giroux, M. L., Fardal, A. A., & Shull, J. M. 1995, ApJ, 451, 477 

Green, J. C., Froning, C. S., Osterman, S. et al. 2012, ApJ, 744, 60 

Gunn, J. & Peterson, B. 1965, ApJ, 142, 1633 

Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Heap, S. R., Williger, G. M., Smette, A. et al. 2000, ApJ, 534, 69 

Hewett, P. G. & Wild, V. 2010, MNRAS, 405, 2302 



- 19 - 

Hodge, P. E. 2011, in Astronomical Data Analysis Software and Systems XX, eds. I. N. 
Evans, A. Accomazzi, D. J. Mink, & A. H. Rots, [A. S. P. Conf. Series 442, ASP, 
San Francisco), 391 

Hogan, C. J., Anderson, S. F. & Rugers, M. H. 1997, AJ, 113, 1495 
Hogg, D. W. 1999, astro-ph/9905116 

Isobe, T., Feigelson, E. D. & Nelson, P. I. 1986, ApJ, 306, 490 

Jakobsen, P., Boksenberg, A., Deharveng, J. M., Greenfield, P., Jedrzejewski, R., & Paresce, 
F. 1994, Nature, 370, 35 

Janknecht, E., Reiniers, D., Lopez, S. & Tytler, D. 2006, A&A, 458, 427 

Kim, T.-S., Carswell, R. F., Cristiani, S., D’Odorico, S., & Giallongo, E. 2002, MNRAS, 
335, 555 

Kim, T.-S., Parti, A. M., Carswell, R. F., & Muller, V. 2013, A&A, 552, 77 
Kirkman, D., Tytler, D., Lubin, D., & Charlton, I. 2007, MNRAS, 376, 1227 
Komatsu, E., Smith, K. M., Dunkley, J. et ah 2011, ApJS, 192, 18 

Kriss, G. A. 1994, in Astronomical Data Analysis Software and Systems III, eds. D. R. 
Crabtree, R. J. Hanisch, & J. Barnes, {A. S. P. Conf Series 61, ASP, San Francisco), 
437 

Kriss, G. A., Shull, J. M., Oegerle, W. et al. 2001, Science, 293, 1112 
Madau, P. & Meiksin, A. 1994, ApJ, 433, L53 
Madau, P. & Rees, M. J. 2000, ApJ, 542, 69 

Martin, D. C., Fanson, J., Schiminovich, D. et al. 2005, ApJ, 619, LI 
McQuinn, M., Lidz, A., Zaldarriaga, M., et al. 2009, ApJ, 694, 842 
McQuinn, M. & Worseck, G. 2014, MNRAS, 440, 2406 
Meiksin, A. 2005, MNRAS, 356, 596 

_ 2006, MNRAS, 365, 807 

_ 2009, Rev. Modern Phys., 841, 1405 (MOO) 





Miralda-Escude, J. 1998, ApJ, 501, 15 

Mortlock, D. J., Warren, S. J., Venemans, B. P. et al. 2011, Nature, 474, 616 
Morton, D. C, 1991, ApJS, 77, 119 

Murdoch, H. S., Hunstead, R. W., Pettini, M. & Blades, J. C. 1986, ApJ, 309, 19 

Osterbrock, D. E. & Ferland G. J. 2006 Astrophysics of Gaseous Nebulae and Active Galactic 
Nuclei (University Science Books, Sausalito, GA) 

Reimers, D., Kohler, S., Wisotzki, L., Groote, D., Rodriguez-Pascual, P. & Wamsteker, W. 
1997, A&A, 327, 890 

Ribaudo, J., Lehner, N. & Howk, J. G. 2011, ApJ, 736, 42 
Richards, G. T., Strauss, M. A., Fan, X. et al. 2006, AJ, 131, 2766 
Schneider, D. P., Richards, G. T., Hall, P. B. et al. 2010, AJ, 139, 2360 
Scholz, T. T. & Walters, H. R. J. 1991, ApJ, 380, 302 
Shapiro, P. R. 1986, PASP, 98, 1014 

Shull, J. M., France, K., Danforth, G. W. et al. 2010, ApJ, 722, 1312 
Shull, J. M., Tumlinson, J. Giroux, M. L. et al. 2004, ApJ, 600, 570 

Smette, A., Heap, S. R., Willinger, G. M., Tripp, T. M., Jenkins, E. B., & Songalia, A. 
2002, ApJ, 564, 542 

Steinhardt, G. L. & Elvis, M. 2011, MNRAS, 410, 201 

Stengler-Larrea, E. A. Boksenberg, A., Steidel, G. G. et al. 1995, ApJ, 444, 64 
Storrie-Lombardi, L. J., McMahon, R. G., Irwin, M. J., & Hazard, G. 1994, ApJ, 427, L13 
Stromgren, B. 1939, ApJ, 89, 526 

Syphers, D. 2010, PhD dissertation. University of Washington 
Syphers, D., Anderson, S. F., Zheng, W. et al. 2009a, ApJS, 185, 20 

_ 2009b, ApJ, 690, 1181 

_2011a, ApJ, 726, 111 



_2011b, ApJ, 742, 99 

_2012, AJ, 143, 100 

Syphers, D. & Shull, J. M. 2013, ApJ, 765, 119 
_2014, ApJ, 784, 42 

Teller, R., Zheng, W., Kriss, G. A., & Davidsen, A. F. 2002, ApJ, 565, 733 
Tittley, E. R., & Meiksin, A. 2007, MNRAS, 380, 1369 
Trainor, R. F. & Steidel, C. C. 2012, ApJ, 752, 39 
Tytler, D. 1987a, ApJ, 321, 49 
_ 1987b, ApJ, 321, 69 

Vanden Berk, D. E., Richards, G. T., Bauer, A. et al. 2001, AJ, 122, 549 

Veron-Getty, M.-P. & Veron, P. 2010, A&A, 518A, 10 

Venemans, B. P., Bahados, E., Decarli, R. et al. 2015, arXiv 15020.1927 

Venemans, B. P., Findlay, J. R., Sutherland, W. J. et al. 2013, ApJ, 779, 24 

Vestergaard, M., & Peterson, B. 2006, ApJ, 641, 689 

White, R. L., Becker, R., Fan, X., & Strauss, M. A. 2003, AJ, 126, 1 

Wilson, J. G., Henderson, G. P., Herter, T. L. et al. 2004, in Ground-based Instrumentation 
for Astronomy, Proc. SPIE 5492, eds. A. F. M. Moorwood & M. lye, 1295 

Worseck, G., Prochaska, J. X., Hennawi, J. F. & McQuinn, M. 2014, arXiv 1405.7405 

Worseck, G., Prochaska, J. X., McQuinn, M. et al. 2011, ApJ, 733, 24 

Wu, X., Wang, F., Fan, X. et al. 2015, Nature, 518, 512 

York, D. G., Adelman, J., Anderson, J. E. Jr. et al. 2000, AJ, 120, 1579 

Zheng, W. & Davidsen, A. F. 1995, ApJ, 440, L53 

Zheng, W., Kriss, G. A., Deharveng, J.-M. et al. 2004, ApJ, 605, 631 

Zheng, W., Kriss, G. A., Teller, R. G., Grimes, J. P. & Davidsen, A. D. 1997, ApJ, 475, 469 



- 22 - 

Zheng, W., Meiksin, A., Pifko, K. et al. 2008, ApJ, 686, 195 


This preprint was prepared with the AAS lATpX macros v5.2. 



23 


APPENDIX A 
ACS PRSIM SPECTRA 


Approximately half of all the known He II qnasars were discovered in three ACS prism 
snrveys (GO 10907, 11215, 11982: PI Anderson; Svphers et ah 2009a : Svphers et ah 2009b). 
The ACS PR130L prism offers a low spectral resolntion, R 170 — 380 between 1360 and 
1250 A. While the potential merit of these prism spectra is limited, their sample size is 
signihcant. It is therefore interesting to explore the proximity prohles in the “other half” of 
the He II qnasar sample. We estimated the proximity prohles in the pnblished prism spectra 
of 15 qnasars at z < 3.6. The spectra of six other qnasars at higher redshifts were not nsed 
becanse the spectral resolntion degrades signihcantly redward of 1400 A. 


To nnderstand the effect of spectral deconvolntion, we simnlated thonsands of spectra 
with different proximity prohles (0 — 25 A) at 2 ; = 3.5 and a pixel scale of 0.1 A, then 
convolved with a Ganssian kernel appropriate to the wavelength-dependent ACS resolntion. 
The inpnt spectra were of simple constant hnx, with a proximity prohle that starts at 
1367 A and declines linearly towards zero at the end of the proximity zone. The hnx at 
shorter wavelengths of the proximity zone was set to zero. The noise level was set for an 
exposnre time of 4500 sec, which is typical for these prism observations. As the S/N level 
in prism spectra is qnite high (> 10), onr simnlation resnlts are not sensitive to noise. Each 
of these simnlated spectra was smoothed to a resolntion of the PR130L prism and then 
binned with small wavelength ohsets. These ohsets rehect the possibility that a spectrnm 
may shift within a prism-data pixel and snbseqnently ahect the deconvolntion resnlt. These 
simulated prism data were deconvolved using the IDL task MAX_LIKELIH00D with a point- 
spread-function kernel of two pixels and then compared with the input data. 

Our measurements of proximity zones in approximately 90% of the simulated spectra 
were within an error of one prism pixel (~ 4 A), and the others recovered within two pixels. 
A conservative estimate for the errors in prism data as 7 A (~ 5 Mpc at z = 3.5), which are 
considerably higher. 

We compared the COS/G140L spectra of quasars SDSS1253-I-6817 and SDSS2346—0016 
with their counterparts in the ACS/PR130L prism data. As shown in Figure El prism spectra 
display an extended wing at a level < 10% of the unattenuated flux. This may be due to 
the extended wing of the instrumental line-spread function or potential misalignment of a 
few individual prism images. We measured the proximity zone to a flux level of 10% level. 
The proximity prohles in the prism data are consistent with that in the COS/G140L data: 
SDSS1253-I-6817 displays a broad prohle and SDSS2346—0016 does not. 

Most of the ACS prism data yielded no detection of proximity prohles (< 5 Mpc). In 
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Table H] we list three quasars whose proximity prohles in the ACS prism spectra are larger 
than 7 A. They are at relatively high redshifts, and the data suggest that the decline of 
proximity-zone size may start at 2 ; ~ 3.5. 

APPENDIX B 

PROXIMITY EFFECT in SDSS DATA 


The study of the proximity effect has been based on high-resolution quasar spectra 


(Murdoch et al. 

1986; 

Baitlik. Duncan & Ostriker 

1988; 

Dall’Adio. Wisotzki. & Worseck 


2008bl) . where the numbers of forest lines are compared in the quasar’s vicinity and at large 
distances. The effect o f enhanced radiation in a quasar’ s vicinit y may be seen even at a 
low spectra resolution. iDaH’Aglio. Wisotzki. fc Worseck fj2008al) used a flux transmission 
technique on the VLT/FORS2 spectra {R ~ 800) of 17 bright quasars at 2.7 < z < 4.1 to 
measure the effective optical depth along the lines of sight and detect a proximity effect. 


The effective optical depth accumulated from most IGM components is estimated as 

r = 0.0021(1 + ;2)3-7 (4) 


flMeiksin 1120061: iKirkman et al.ll2007l: iBecker et al.ll2013[) . Approximately half of this term is 
from those at column density N > 10^^ cm“^, most of which cases can be identihed at the 
SDSS spectral resolution and S/N fl l2.4p . From a power-law distribution of dn/ dNnr oc 
where N is the colum n density and n the number of individual absorbers flTvtlerl Il987a : 


Janknecht et al.l 120061) . we estimated approximately six strong absorbers over a range of 10 


Mpc 66 A on the SDSS wavelength scale). These lines affect roughly 60% of the pixels 
in this bin, assuming that each line affects hve SDSS spectral pixels. The fluxes in other 
pixels in this bin are believed to bear the signature of weak IGM components. Since these 
lines are unsaturated, their intensity is sensitive to the proximity effect. By measuring the 
flux decrements in these pixels, we tried to detect the signature of proximity zones. 


We generated a set of absorption templates r = 0.0021(1 -|- /{I -|- (Rp/rY), where 

r is the distance to the quasar and Rp the proximity size. The range of Rp is 5—35 Mpc. 
The htting task specfit allows such a “userabs” component in the form of opacity vs. 
wavelength. We ran specfit on SDSS spectra with the following components: an underlying 
power law, Gaussian emission components of Ly a (narrow -|- broad), N V, O I, Si IV, G IV, 
He II, and a set of absorption templates. The htting windows consists of two parts at the 
following restframe wavelengths: (1) 1216—1680 A, which cover the red wing of Ly a and 
hve other emission lines fl l2.ip and (2) a number of small windows between 1200 and 1216 
A, where no strong forest line is detected. 


We ran tests on the SDSS spectra two quasars with diherent He II proximity zones. For 
SDSS1253-I-6817, our htting results suggest a symmetrical Ly a prohle, implying a strong 




































proximity profile. The average opacity Te// = 0.03±0.03 over a bin 20 —80 A blueward of the 
Ly a wavelength. The calculation includes 22 pixels that are not in absorption windows. For 
SDSSJ1711+6052, the htting result suggests that the blue Ly a wing is weaker than the red 
wing, as expected from IGM absorption. The average optical depth in a similar wavelength 
range is Te// = 0.46 ± 0.33. The S/N in the second quasar is ~ 7 per pixel in the proximity 
zone, and the errors in our measurements are dominated by pixel-to-pixel variations. A 
higher average optical depth with noticeable pixel-to-pixel fluctuations is consistent with a 
signihcant structure in the forest lines and a weak proximity effect. These fitting results 
suggest that the signature of a broad proximity zone may be detected in SDSS spectra. 



Table 1. Summary of Observations 


Quasar 

R.A. 

Decl. 

Observation Exposure Time 

Grating Gentral 


(J2000) 


(sec) 

Wavelength (A) 


Data of GO 12249 

SDSS1253+6817 

12 53 53.715 

68 17 14.20 

2011 May 

14095 

1105/1280 

SDSS2346-0016 

23 46 25.662 

-00 16 00.47 

2010 Nov., Dec. 

20737 

1105/1280 

SDSS1711+6052 

17 11 34.412 

60 52 40.39 

2011 Apr., Oct. 

23950 

1105 

SDSS1319+5202 

13 19 14.205 

52 02 00.11 

2011 May 

26643 

1105 

Archival Data of GO 11742, 13013 

SDSS1237+0126 

12 37 48.993 

01 26 06.90 

2010 Jun. 

6212 

1105 

QSO2149-0859 

21 49 27.770 

-08 59 03.61 

2013 Apr. 

7561 

1105 

SDSS0936+2927 

09 36 43.511 

29 27 13.60 

2011 Jan. 

4739 

1105 

QS01630+0435 

16 30 56.340 

04 35 59.42 

2013 Apr. 

7908 

1105 

QSO2157+2330 

21 57 43.630 

23 30 37.34 

2013 Jul. 

8074 

1105 

QS00233-0149 

02 33 06.010 

-01 49 50.58 

2013 Aug. 

2021 

1105 

QS00916+2405 

21 57 43.630 

23 30 37.34 

2013 Dec. 

3075 

1105 



Table 2. Proximity Zone Measurement 


Quasar 

Redshift' 

This work 

a 

SDSS 

Observed Gontinuum Flux /ao(304) 
lO””” erg s“’ cm“^ A 

Magnitude 

A7i45o 

Proximity Zone Size 
A Mpc 


Data of GO 12249 

SDSS1253+6817'’ 

3.476 

±0.004 

3.4727 

19 

-27.3 

20 

12 

SDSS2346-0016" 

3.511 

±0.003 

3.4895 

19 

-28.0 

6 

3.5 

SDSS1711+6052^ 

3.823 

±0.007 

3.8269 

8 

-26.6 

5 

2.5 

SDSS1319+5202^ 

3.93 

±0.01 

3.8991 

2 

-28.2 

10 

4.7 

Archival Data of GO 11742, 13013 

80880936+2927” 

2.9253 

; ± 0.004 

2.9239 

15 

-27.4 

20 

16 

QSO2157+2330 


3.142'’ 


15 

-27.7 

26 

19 

SDSS1237+0126'’ 


3.1448'’ 


15 

-26.8 

11 

8 

QSO2149-0859 


3.259'= 


6 

-26.9 

4 

2.7 

QS00233-0149 


3.314'= 


11 

-27.2 

11 

7.2 

QS00916+2405 


3.440'= 


18 

-27.0 

7 

4 

QS01630+0435 


3.788'= 


23 

-28.3 

10 

5.1 


“^The nominal redsliift uncertainty is 0.002 for SDSS spectra. 

’’Optical spectra checked. Strong H I absorption found near the end of the proximity zone. 
'’Optical spectra checked. No strong H I absorption found near the end of the proximity zone. 
'’O I A1304 emission is contaminated by absorption. The redshift value is from 
®http://www.stsci.edu/hst/phase2-public/13013.pro 


Worseck et ah (^20141 ). 
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Table 3. He II Measurement Based on Published COS and STIS Data 


Quasar RedshifP^ Magnitude Proximity Zone Size Reference 

Tfl450 (^) (Mpc) 


HS1700+6416'^ 

2.751 ±0.003 

-29.4 

4 

3.7 

Q1216+1656 

2.818 

-27.3 

15 

13 

HS1024+1849'" 

2.8475 

-27.0 

21 

18 

4C57.27 

2.858 

-27.1 

11 

9 

HE2347-4342" 

2.885 ±0.005 

-29.3 

0 

0 

SDSS1508+1654'^ 

3.1716 

-27.3 

18 

12 

PKS1935-692'’ 

3.185 

-28.3 

25 

17 

SDSS0856+1234'^ 

3.1948 

-27.0 

20 

14 

SDSS0955+432 ^ 

3.2388 

-26.2 

18 

12 

QS00302-003^ 

3.2860 

-28.3 

17 

11 

SDSS0915+4756" 

3.3369 

-27.9 

6 

3.9 

SDSS2345+1512'^ 

3.5880 

-26.6 

5 

2.8 

SDSS2257+0016^ 

3.7721 

-27.2 

5 

2.6 


Svnhers fc Shull 120131 
Svphers et al. 120121 
Svnhers et al. 120121 
Svphers et al. 120121 
Shull et al. 120101 
Svphers et al. 120121 
Anderson et al. 119991 
Svphers et al. 120121 
Svphers et al. 120121 
Svphers fc Shull 120141 
Svnhers et al. 120121 
Svnhers et al. 120121 
Svnhers et al. 120121 


^From the SPSS D R7 quasar cat alog (Schneider et ah 2010ll . except: HS1700+6416 


llTrainor fc Steidel 112012); HE2347—4342 flReimers et al 


2014 1: PKS1935-692 flAnderson et al.l 1199911 : 


1199711 : QS00302 —003 llSvphers fc Shull 
Q1216+1656, 4C57.27 f Veron-Cettv fc Veron 


201011 ■ The nominal redshift uncertainty is 0.002 for SDSS quasars. 


'^Optical spectra checked. Strong H I absorption found near the end of the proximity zone. For 
HS1700+6416, this line is present at z = 2.744. 


Optical spectra checked. No strong H I absorption found near the end of the proximity zone. 
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Table 4. He II Measurement Based on Published ACS Prism Data 


Quasar 

Redshift^ 

Magnitude 

d4i45o 

Proximity Zone Size 
(A) (Mpc) 

Reference 

SDSS1042+5129 

3.3864 

-27.1 

26 

17 

Svnhers et al. f2011b') 

SDSS1007+4723 

3.4084 

-26.2 

9 

6 

Svnhers et al. f2009b') 

SDSS1442+0920 

3.5286 

-28.3 

13 

8 

Svnhers et al. f2009b') 


'^From the SDSS DR7 quasar catalog ( Schneider et al. 20101) . 
certainty is 0.002 


The nominal redshift un- 
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Figure Captions 


Fig. 1.— HST/COS G140L spectra of four quasars. The data are binned by six pixels 
0.5 A, 0.85 resolution elements). The green dotted curves are the errors, and the red 
dashed lines mark the wavelengths of He II Ly a. Vertical bars mark the wavelength of 
strong absorption lines identihed from the SDSS spectrum counterparts, and arrows mark 
the edge of proximity zones. 

Fig. 2.— Near-infrared spectrum of SDSS2346—0016, taken with the ARC 3.5-m telescope 
and TripleSpec instrument. Gaps in wavelengths are due to the removal of poor data in 
the bands of severe atmospheric absorption. The inset displays the wavelength region for 
redshifted Mg II emission. We used it to determine the systemic redshift oi z = 3.511. 

Fig. 3.— HST/COS G140L spectra of seven other quasars in the HST archive, plotted 
in the quasar’s restframe. The data are binned by six pixels (~ 0.5 A). The green dotted 
curves are the errors, and the red dashed lines mark the wavelengths of He II Ly a. The 
Earth symbols mark the wavelengths of the strong geocoronal lines around 1216 and 1302 A, 
and a weak one at 1356 A. Vertical bars mark the wavelength of the strong absorption lines 
identihed from the SDSS spectrum counterparts, and arrows mark the edge of proximity 
zones. 


Fig. 4.— Normalized spectra of quasar SDSS1253-I-6817. The lower panel is the He II Ly a 
region in the observer’s wavelength frame (COS/G140L), and the upper panel is the aligned 
H I counterparts (SDSS). The proper distance scale in the top label is appropriate for both 
panels. The green shaded region marks the wavelength range for the He II Ly a break, given 
the redshift uncertainty {z = 3.472 — 3.480). The magenta arrow marks the wavelength for 
the ending point of the estimated proximity zone. The blue lines mark respective baselines 
and data errors. The red bars mark the strong absorption features identihed in the SDSS 
spectrum. 


Fig. 5.— Size of proximity zone vs. redshift. In the upper panel, red boxes are COS 
data, GO 12249; magenta boxes: archival data of GO 11742 and 13013; green boxes: COS 
data from literature; the cyan star: STIS data of PKS1935—692. The downward arrow is 
for la limit with n on-de tection. The lo wer panel is for H I. Black cro sses are from Table 


1 of Carilli et a 


Mortlock et ah 


fl2010h: green ci r cles: IVenemans et ah fj2013l. l2015l) : and red triangles: 


(120111) ; IWu et ah (120151) . The redshift ranges in the two panels are scaled 


to represent the same ratio (1.0:0.65) of cosmic ages from the Big Bang. 


Fig. 6.— Size of proximity zone vs. absolute magnitude. The symbols are the same as Figure 
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ls] Yellow curves represent model predictions of Rg oc The very luminous quasars 

(M < —29) show proximity-zone sizes that are considerably smaller than anticipated. 

Fig. 7.— Simulated He II ionization structure, assuming a homogeneous IGM. The lower 
panel; log(He II fraction), middle panel; gas temperature in units of 10^ K, and upper 
panel; normalized flux. The quasar He II ionization rate has been boosted to 2 x 10^^sec“^. 
This extreme case assumes the quasar has been observed in a relatively quiescent state. 
The ionization levels assume the boosted rate, as it would take about 9 Myr for ionization 
equilibrium to be established at the lower rate at the ionization front. Five curves represent 
an expansion time sequence of 10 (cyan), 20 (blue), 30 (green), 50 (magenta), and 100 
(red) Myr after the quasar turns on. The age of the quasar producing the ionized zone will 
generally be smaller (see text). The proximity prohle in SDSS1253-I-6817 is overplotted. 

Fig. 8.— Comparison of ACS/PR130L and COS/G140L spectra of two quasars. In the top 
panels, the green histograms are the HST data, and black histograms the deconvolved data. 
The proximity zones are marked between the respective magenta arrow and red vertical line, 
which is at the wavelength of He II lya. The zero-flux baselines are marked by blue dashed 
lines. 
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Fig. 5. 


















-37- 



-29 -23 -27 -26 -25 


Magnitude 


Fig. 6. 












log Hell Fraction Temperature (10^) Flux 


38 



5 10 15 20 25 


Proper Distance (Mpc) 
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